Dynamic crossover in the global persistence at criticality 
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We investigate the global persistence properties of critical systems relaxing from an initial state 
with non- vanishing value of the order parameter (e.g., the magnetization in the Ising model). The 
persistence probability of the global order parameter displays two consecutive regimes in which 
it decays algebraically in time with two distinct universal exponents. The associated crossover is 
controlled by the initial value mo of the order parameter and the typical time at which it occurs 
diverges as mo vanishes. Monte-Carlo simulations of the two-dimensional Ising model with Glauber 
dynamics display clearly this crossover. The measured exponent of the ultimate algebraic decay is 
in rather good agreement with our theoretical predictions for the Ising universality class. 
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INTRODUCTION 



In spite of many efforts during the last decades, a de- 
tailed description of the non-equilibrium dynamics of sta- 
tistical systems is still lacking. In addition to disordered 
and glassy systems [l[ even simpler ones (e.g., pure mag- 
nets) have revealed unexpected non-equilibrium dynami- 
cal behaviors such as aging 0,0]. Close to critical points, 
some aspects of these collective behaviors become largely 
independent of the microscopic details of the system (uni- 
versality) and therefore they can be very effectively inves- 
tigated by means of simplified models, among which field- 
theoretical (FT) [| ones. Within this approach one stud- 
ies, for instance, the non-equilibrium relaxation from a 
state with non- vanishing value too of the time-dependent 
average order parameter m(t) (e.g., the magnetization 
of a ferromagnet) . It turns out [4[ that, after a non- 
universal transient, m(t) grows in time as m(t) oc mot 
for t <C T m oc m Q 1 ' k whereas, for t » r m , m(i) decays 
algebraically to zero as to (t) oc t-PIW. These differ- 
ent time dependences are characterized by the univer- 
sal exponents 9' (the so-called initial-slip exponent Q) 
and k = 9' + [3/(vz), where f3, v and z are the usual 
static and dynamic (equilibrium) critical exponents, re- 
spectively. In addition to the actual time dependence of 
m(t), the 'persistence properties of thermal fluctuations 
5m(t) around m(t) provide useful informations on the 
dynamics of the system (especially on its non-Markovian 
nature) and indeed in recent years they have attracted, 
in various contexts, considerable attention both theoret- 
ically H and experimentally [(|. For a stochastic pro- 
cess {X(t)}t>o with zero mean, the persistence proba- 
bility P(t) is defined as the probability that X{t) does 
not change sign up to time t. A related quantity which 
might be used to characterize the process is the first- 
crossing (or first-passage) time t CIO , defined as the time 



it takes X(t) to change its sign for the first time. Clearly, 
the probability distribution p(t cm ) is related to P(t) via 
p(t C ro) = —P'(t cro ). For a ferromagnet at finite temper- 
ature T the persistence probability of a single spin (i.e., 
of the fluctuations of the local order parameter) decays 
always exponentially in time due to fast thermal fluctua- 
tions. The long-time decay of the persistence probability 
P c {t) of 8m(t) (i.e., of the fluctuations of the global order 
parameter), instead, becomes algebraic upon approach- 
ing the critical point: P e (t) oc t~ e , with a non-trivial ex- 
ponent 9 > 0. This exponent was calculated analytically 
for some ferromagnetic models with different dynamics 
(Models A and C 0|) and random initial condition [1,0] 
(i.e., mo = 0) and for reaction-diffusion systems in the di- 
rected percolation (DP) universality class, initially in the 
active phase [10( (i.e., mo = 1 on the lattice, where the or- 
der parameter is the spatial density of reacting particles). 
We shall denote 9q = 9(mo — 0) and 6*00 = 9(mo — 1). In 
the previous cases analytical progresses rely on the fact 
that the global fluctuating order parameter m + 5m has 
a Gaussian distribution for all finite times t in the ther- 
modynamic limit. Indeed, for a d-dimensional system of 
linear size L, to + 5m is given by the sum of L d random 
variables (the local fluctuating degrees of freedom, e.g., 
spins in ferromagnets) which are correlated only across 
a finite, time-dependent correlation length In the 

thermodynamic limit L/£(t) 3> 1 the central limit theo- 
rem implies that to + <5to is a Gaussian process, for which 
powerful tools have been developed in order to determine 
the persistence exponent [B|, [ll[ . Remarkably, under the 
additional assumption that the process is Markovian, 9 
can be related to known critical exponents, via two (hy- 
per)scaling relations 6*o = no = —9'+(l—r}/2)/z (where rj 
is the Fisher exponent) and 6^ = fi^ = l + rf/(2z) which 
are valid for spin systems and for the DP universality 
class, respectively. For toq = (weak) non-Markovian 



2 



corrections were found at two-loop and one-loop order 
in a dimensional expansion around d = 4, for Model A 
and Model C, respectively, in rather good agreement with 
Monte Carlo estimates Oo = 0.237(3) in dimension d = 2 
Qjl and O ^ 0.41 in d = 3 (l3f. For completely ordered 
initial conditions (mo = 1), instead, these corrections 
have not been investigated so far beyond the case of DP, 
where they appear at one-loop order [lol ] and they were 
measured numerically in d = 1 [3| • These results imply 
that the average and variance of the first-crossing time 
t CTO are both finite for DP in d < 4 but not for the spin 
models studied in Refs. [9l. HlT|. 

The evidence accumulated so far in two different mod- 
els with mo = and 1 show that the critical persistence 
properties actually depend on mo. It is therefore natu- 
ral to investigate how the corresponding crossover occurs 
when mo is varied within the same model [3| and which 
arc the associated universal features. 



m is 



ip(x,t) = <p(x,t)-M(t) , where M(t) = {(f(x,t)) (3) 

is the average magnetization and (...) stands for the av- 
erage over the possible realizations of the stochastic noise 
£. Defining M(t) as 



M(t) = j- d ld d x^( Xl t) 



(4) 



we are interested in the probability P c (t) that M does 
not change sign in the time interval [0, t] following the 
quench to T c . 

We first present the result of the Gaussian approxi- 
mation which is exact in dimension d > 4 and which 
is obtained by neglecting non-linear terms in ip in the 
Langevin Eq. {T]) expressed in terms of ip(x, t) and m 2 = 
g M 2 /2 (see, e.g., Ref. Q): 



UNIVERSAL SCALING BEHAVIOR OF THE 
PERSISTENCE PROBABILITY 

Here we focus primarily on the Ising model (on a d- 
dimensional hypercubic lattice) with Glauber dynamics 
quenched to its critical point and we study its persis- 
tence properties both analytically and numerically. The 
universal aspects of the relaxation of this model are cap- 
tured by the so-called Model A [7J for the n-component 
fluctuating local order parameter (f(x, t) (e.g., the coarse- 
grained density of spins in the Ising model): 



■nd t <p(x,t) = - 



SH[(p] 
6ip(x,t) 



C(x,t) 



(1) 



where ((x, t) is a Gaussian white noise with (((x, t)) = 
and (C(x,t)C(x',t')) = 2r]T6(x - x')5{t - /), r\ is the 
friction coefficient (set to 1 in the following) and T the 
temperature of the thermal bath. In Eq. |TJ, 7i is the 
0(n)-symmetric Landau-Ginzburg functional which, for 
n = 1, describes the universal aspects of the static prop- 
erties of the Ising model: 



H[p] = J d d x 



1 2 

2 r0lf 



> 2 ) 2 



(2) 



where ro is a parameter which has to be tuned to a critical 
value ro jC when approaching the critical temperature T c 
(here ro, c = 0) and go > is the bare coupling constant. 

At the initial time t = 0, the system is in a random 
configuration with mean magnetization [ip(x, 0)]j = Mq 
and short-range correlations [<p(x, Q)ip(x', 0)], = t ~ 1 S(x— 
x') ([. . .]i stands for the average over the distribution 
of the initial configuration), where Tq 1 is irrelevant in 
determining the leading scaling properties [H (therefore 
we set Tq -1 = 0). For n = 1 (the case n > 1 is shortly 
discussed below) , the stochastic process we are interested 



d t ip{x,t) = [V 2 -m 2 (t)]ip(x,t) +((x,t) 

1 (5) 
where d t m(t) + -m (t) — . 

The equation of motion for M [see Eq. (jJJ)] readily 
follows: d t M(t) + m 2 {t)M{t) = C(i), where C(t) = 
L~ d Jd d x ((x, t) is still Gaussian and, from Eq. |5]), 
m 2 (t) = (2i/3 + Wq 2 )" 1 where m(t — 0) = mo. In 
the limit t <C m^ 2 , one finds that M(t) executes a sim- 
ple random walk and therefore P c (t) oc t~ 1 / 2 [3]. In 
the opposite limit t 3> Tig -2 , M(t) satisfies the Langevin 
equation d t M{t) + (3/2t)M(t) = ((t) which turns into 
d T M{r) = t){t) when the variables M{t) = t 3/2 M(t), 
T](t) = i" 3/2 C(i)/4 and r = t 4 are introduced 0. Ac- 
cordingly, P c (t) cx r -1 / 2 , i.e., P c (t) oc i~ 2 , which de- 
cays to zero more rapidly than for t "i^ 2 and clearly 
shows the existence of two different regimes. The full 
crossover function can be obtained by noticing that M(t) 
is a Gaussian process and therefore it is completely de- 
termined by its two-time correlation function CjrAt, t') — 
(M(t)M(f)) which can be easily calculated 0. Fol- 
lowing Ref. Q, one introduces the normalized process 
X(t) = Af(t)/(M 2 (i)) 1 / 2 , the two-time correlation func- 
tion of which is given by (i > t') [lc| : 



Lit') 



^/C^(t,t)C^(t',t') L(t) (6) 



(X(t)X(t')} 



where L{t) = yj (t + l) 4 - 1 

and the dimcnsionlcss time variables are given by t — 
t/r m with t,„ = 3/(2mo). In terms of the logarithmic 
time T = In L(t), X(T) is a stationary Gaussian pro- 
cess with purely exponential correlations (X(T)X(T')) — 
exp (— \T — T'|), for which the persistence probability is 
known exactly p, [l3|- In the limit t,r m ^> £ m icr (where 
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^micr is some non-universal microscopic time scale) one 
finds 



Pc{t) OC 



L\txnicT ) _ {j~m / ^micr) 



-1/2 



L(t) y/(l + t/r m y - 1 

(V* , t < T m 
It" 2 , t > T m 



(7) 



in agreement with the results obtained by mapping the 
process on a random walk. The mean first-crossing time 
£ CT o resulting from P c (t) is finite for too 7^ 0, whereas 
its variance is not. Although the Gaussian approxima- 
tion is exact only for d > 4, it clearly displays for fi- 
nite mo (i.e., finite r m ) an interesting crossover - 
which is expected to occur for generic d above the lower 
critical dimension of the model — between two succes- 
sive regimes in which P c (t) decays algebraically with two 
different exponents, #0 and 001 respectively. The first 
regime (t <C r TO ) has a temporal extent ~ r m which 
increases upon decreasing too and indeed, in the limit 
mo - * considered in Refs. [il, Q , it is the only one ac- 
cessible. The crossover to the second regime takes place 
for t ~ T m , with T m decreasing upon increasing mo 7^ 
and indeed r m tends to zero in the case too — ► 00 consid- 
ered in Ref. for the DP. 

To go beyond the Gaussian approximation, we take 
advantage of recent results about the aging behavior of 
systems described by Eqs. (JlJ and {2J following a quench 
from a state with mo ^ to the critical point [l5j |. 
As mentioned in the introduction, M(t) (and therefore 
X(t) = M(t)/(M 2 (t)y/ 2 ) is a Gaussian process in the 
thermodynamic limit. Taking into account the scaling 
form of Cjfr(t,t') discussed in Ref. [la], one expects for 
t>t': 



(X(t)X(t')) = (t/t')-^F x (t'/t,t/T m ) 



(8) 



where r m = (i3 m TOo) _1 ^ K , B m is a non-universal constant 
which can be fixed according to some normalization con- 
dition (see Ref. [IH for details), Fx is a universal scaling 
function (with Fx (1, J/) = 1) and ^0 = -0' + U-rj/2)/z. 
For t/r m <C 1 one recovers the result of Refs. [alSl- How- 
ever, a non- vanishing mean too of the initial order param- 
eter affects the behavior of the temporal correlations as 
soon as t ~ r rn and in particular for t, t' 3> r m . The 
limiting behaviors of Eq. ([8]) can be derived from those 
of Cjfr(t, t') discussed in Ref. [HI]: 



(X(t)X(t')) 



(t/t')-^f Q (t'/t), t'<t«T„ 
{t/t')-^ foo{t' /t) , i>i'>T„ 



(9) 



where f Q (x) = Fx^O), /^(a;) = x"»-M- J F A .( iB> y -» 
00) with fo.no (0 ) fi nite [/o,oo(l) = 1] and, for d < 4, 
A*oo = 1 + d/(2z)[25j. Note that if /o, 00 (x) is constant for 
x £ [0, 1], (X(t)X(t')) can be cast asymptotically in the 
form © with a suitable choice of L(t) and therefore the 



process X(t) can be mapped into a Markovian one with 
Pc{t micl - < t < r m ) - and P c (t > r m ) - t~^. In 
the scaling limit, the relevant time scales for the problem 
of persistence are the time t elapsed since the quench and 
the (mesoscopic) time scale r m set by the value of the ini- 
tial magnetization. Accordingly, one expects (also on the 
basis of the Gaussian approximation) a scaling behavior 



p c (t) = A P r 9o v{t/T m ) , v(x > 1) 



(10) 



where Ap is a non-universal constant such that V(0) = 1, 
V(x) is a universal scaling function, r m is fixed by the 
behavior of the magnetization as in Ref. [151 ]. and #00 
is a yet undetermined exponent. Equation (|10|) implies 



(£ C ro) = Jo dtP c (£) ~ T m ~ e ° and therefore the first- 
crossing time i cro has a finite average for #0 < 1 < $00 and 
a finite second moment (^ ro ) = 2 J °°dttP c (t) ~ r^j~ e ° for 
6>o < 2 < 6*00. Note that #0.00 = Mo. 00 for a Markovian 
process. The corrections to #0 due to non-Markovian 
effects (reflected in the fact that fo y0 o(x) in Eq. ((9]) are 
not constant) were obtained at two-loop order in Ref. 0] . 
For these corrections have not been calculated. They 
can be determined by focusing on the case too — > 00 of 
Eq. (O, i.e., on foo(x), which was calculated in Ref. [HI 
up to one loop in the e-expansion (e = 4 — d). Taking 
advantage of the results therein one finds, in logarithmic 
time T = lni (T > T'), 

(X(T)X(T')) = c-^ T - T 'U(e-( T - T '>), (11) 
where A{x) = 1 + eAi(x) + 0(e 2 ) , 

Ai(x) = 2[fc(x) - /c(l)]/3 and f c (x) is given in 
Eq. (B.31) of Ref. [H]. The universal function Ai{x) 
carries the signature of the non-Markovian corrections to 
#00 which, in turn, can be calculated by using the per- 
turbation theory of Ref. [9j (see Eq. (14) therein): 



n = 



= 1-6 



2fi 



= l + e0.0273.. + O(e 



, x H 
ax- — 

la (1 

2\ 



.)3/2 



0{e 2 ) 
(12) 



where ^ = l + d/(2z) = 2(1 -e/8) + 0(e 2 ) and the inte- 
gral has been computed numerically. At variance with 9q 
for Model A, non-Markovian corrections appear alr eady 
at one-loop order, as in the DP and Model C fiof . 
Compared to its Markovian and Gaussian approxima- 
tions, the resulting exponent 6>oo = 2 — e 0.195.. + 0(e 2 ) 
is respectively increased and decreased. Accordingly, t CYO 
has a finite average (t cro ) ~ tIq ^ 8o '/ K which diverges 
with a universal exponent upon reducing too. In the op- 
posite limit too — ► 00, instead, (tcro) is finite and deter- 
mined by imicr (as for the DP }10|). 

It is also instructive, though less relevant for the de- 
scription of systems other than lattice models, to calcu- 
late the persistence probability of a n-component vec- 
tor order parameter (p = {(p±, (f2, .., tp n ) with n > 1 
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and Model A dynamics (see Eqs. ([1} and J2])). When 
the initial magnetization mo is finite, fluctuations which 
are parallel tp <T (x > t) and transverse ip n (x,t) (with 0(n — 
1) symmetry) to the average order parameter m(t) — 
(tp(x,t)) are expected to have distinct persistence prob- 
abilities, P° (t) and Pc(t), respectively. Taking advan- 
tage of the results of Ref. [18( one finds that both of 
them exhibit a crossover as in Eq. (|10|l with n-dependent 
exponents #o( n ) = ^o( n ) (f° r * ^ r m the O(n)- 
symmetry is restored and the results of Ref. [9J apply) 
and (n) which in the Markovian approximation read 
Hi, = 1 + d/(2z) and /i^ = 1 - d/(2z) + 2{3/{vz). Non- 
Markovian corrections contribute already at one loop, 
yielding Vf = 1 + e(0. 115.. + n 0.131.. )/(8 + n) + 0{e 2 ) 
and W = 1 + e 0.055. ./(8 + n) + 0(e 2 ) (an expansion of 
TZ* around d = 2 can be obtained from the analysis of 
Ref. (HI). Interestingly enough, by taking advantage of 
the results of Ref. [18| . P£(t) can be computed exactly 
for n — > oo and 2 < d < 4. In the short-time regime 
t -C r m , one recovers #o( n ~~ * 00 ) = — 2)/4, i.e., the 
same expression as for the case mo = Q whereas, for 
t ^> T m , one finds #£o(ft — > oo) = d/4. In both cases there 
are no non-Markovian corrections and indeed the correla- 
tion function of the normalized process X™ (t) associated 
to M r (t) = L~ d $A d xi>K{x,t) takes the form © with 
L(t) = t( £i - 2 )/ 4 [t + d/(d-2)] 1 /2 an d Tm = 3(d-2)/(4m§), 
allowing also the calculation of the full cross-over func- 
tion. In passing we mention that a similar crossover 
occurs also in the spherical model with purely dissipa- 
tive dynamics: By taking advantage of the results of 
Ref. [20j we find 9oc = /ioo = rf/4 + 1 (in agreement 
with the hyperscaling relation) which differs indeed from 
#o = Mo = (d-2)/4| 



MONTE CARLO SIMULATIONS OF THE 
TWO-DIMENSIONAL ISING MODEL 

To test the expected scaling behavior (see Eq. ([TO"])), 
the corresponding crossover as well as our theoretical es- 
timate of 6oo , we performed Monte Carlo simulations of 
the two-dimensional Ising model (with spins Si — ±1) on 
a L x L lattice with periodic boundary conditions and 
Glauber dynamics. The system is prepared in a random 
configuration with N + up and AT_ down spins, where 
N± = L 2 (l ± mo)/2. At each time step, one site is ran- 
domly chosen and the move Si i— > — s, is accepted or 
rejected according to Metropolis rates. One time unit 
corresponds to L 2 such time steps. The computation of 
the persistence probability P c (t) requires the knowledge 
of the magnetization M(t) = L~ 2 (^2 i s^, which we aver- 
aged over 2000 thermal realizations. After the quench to 
T c , each system evolves until M(t) = L~ 2 J^i s i ~ M(t) 
first crosses zero. P c (t) is then measured as the frac- 
tion of surviving systems at each time t, computed by 
averaging over 10 5 samples. In Fig. [JJ we plot P c (t) for 
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FIG. 1: Persistence probability P c {i) of the two-dimensional 
critical Ising model with Glauber dynamics, for different ini- 
tial magnetization mo. All the curves decay initially with an 
exponent #o — 0.235 and then, after the crossover, with a 
different exponent 6 X ~ 1.7. 

L = 256 and different values of mo ranging from to 
1. The curves with mo > 0.1 display a clear crossover 
between two different algebraic decays, as anticipated by 
our theoretical analysis. We have carefully checked - 
by computing P c (t) for L ranging between 64 and 512 — 
that the displayed crossover is not affected by finite-size 
effects. For small enough t, P c if) decays with an expo- 
nent = 0.235(5) (fully compatible with the avail- 
able estimate [l2[ ), whereas for larger times the power- 
law decay is faster, with > #0 and 

C C) = 1.7(1) ■ (13) 

An analytical estimate of this exponent can be obtained 
by setting e = 2 in Eq. {TJJ), yielding 6^ loop) ~ 1.61 
which is actually in rather good agreement with 9^ C \ 
As expected from the scaling form (fTO)) , the time r m at 
which the crossover occurs in Fig.[JJincreases as the initial 
magnetization mo decreases and eventually no crossover 
is observed for m = [||. In Fig.[2]we plot P c (t)mQ 6 °^ K 
(oc P c {t)T^) as a function of the rescaled time tm^ K (oc 
t/T m ) in order to highlight the scaling behavior (fTO"]) . The 
exponent k is related to [3 = 1/8, v = 1, 9' = 0.191(3) [3 
and z = 2.1667(5) [2lj via the scaling relation n = 9' + 
f3/{uz) :4], which yields k ~ 0.249. For 9 we use the 
estimate O MC ^ obtained from Fig. [JJ The quite good 
data collapse in Fig. [2] is remarkably obtained without 
fitting parameters and therefore it provides a convincing 
numerical evidence of the scaling behavior (|10|) . 

CONCLUSIONS 

We have studied the global persistence probability 
P c (t) for critical system which are initially prepared in 
a state with short-range correlations and non-vanishing 
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FIG. 2: Plot of P c (t)m Q "^ K as a function of the rescaled 
time tm^ K with 8q — 0.235 and k = 0.249, including the 
data presented in Fig. [T] According to Eq. (|10|l the exponent 
of the ultimate power-law decay is c/oo — 1.7. 



average mo of the order parameter. For different models, 
we have shown that P c (t) exhibits a crossover (appar- 
ently overlooked in previous studies) between two distinct 
power-law behaviors which is described by a universal 
scaling function V (see Eq. (flO)) '). Our Monte Carlo sim- 
ulations of the two-dimensional Ising model with Glauber 
dynamics clearly display this crossover. We have calcu- 
lated non-Markovian corrections to the (universal) per- 
sistence exponent 9^ which characterizes the ultimate 
long-time behavior of P c (t) for mo > and it turns out 
to be in good agreement with our numerical estimates. 
It would be interesting to extend beyond the Gaussian 
approximation the calculation of the scaling function V . 

In this paper we have studied the crossover in the per- 
sistence of the global magnetization M of a d-dimensional 
system at criticality. In addition one might have con- 
sidered the magnetization M su b of a d'-dimensional sub- 
manifold of the original system (0 < d! < d) [22| or 
the magnetization My(t) of a subsystem of finite vol- 
ume V = t d [23|]. It turns out that, depending on d! , the 
long-time decay of the persistence probability of Af s „h 
is exponential, stretched exponential, or algebraic [22], 
whereas the persistence probability of My is algebraic 
for £(i) <C I and exponential for £(i) ^> £. 

The dependence on d' and t of the dynamic crossover 
highlighted in this paper and the interplay among the 
various crossovers surely deserve further investigations. 
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